The future is noisy: on the role of spatial fluctuations in genetic switching 
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A genetic switch may be realised by a certain operator sector on the DNA strand from which 
either genetic code, to the left or to the right of this operator sector, can be transcribed and 
the corresponding information processed. This switch is controlled by messenger molecules, i.e., 
they determine to which side the switch is flipped. Recently, it has been realised that noise plays 
an elementary role in genetic switching, and the effect of number fluctuations of the messenger 
molecules have been explored. Here we argue that the assumption of well-stirredness taken in the 
previous models may not be sufficient to characterise the influence of noise: spatial fluctuations play 
a non-negligible part in cellular genetic switching processes. 

PACS numbers: 87.15.Ya,87.16.-b,05.40.-a,02.50.-r 



The Darwinistic theory of evolution through mutation 
and selection is based on occasional errors in the process 
of biological reproduction |l|. At the same time, a bio- 
logical system has to be locally stable in the sense that 
a mutation occurs only once within a large number of 
generations of a species. Biological systems, faced with 
thermal activation at room temperature therefore have to 
find a way to minimise the influence of noise. Particular 
interest in this concern is focused on genetic switches, a 
relatively simple unit in a biological cell at hand of which 
the effects of noise can be studied. 

Genetic switches constitute a part of the centre of com- 
mand of biological cells. By deciding which of two genetic 
codes within the corresponding DNA section is tran- 
scribed, the switch prompts the production of certain 
molecules inside the cell, and thus it controls the sub- 
sequent reactions and the feedback: the genetic switch 
governs the future state of the cell. In fact, the associ- 
ated interplay of different chemical reactions resembles 
a logical electric circuit, and biochemists have therefore 
coined the notion of genetic circuitry 0, H, M. Genetic 
switches, depending on their functional task, come either 
non-cooperative or cooperative. Roughly speaking, coop- 
erative switches will always be concerned with the control 
of vital processes such as reproduction, non-cooperativity 
is associated with less precision and concerns processes 
like respiration. For a large number of cellular systems, 
the steric and kinetic aspects of the biochemistry of the 
molecules involved in genetic switching have been ex- 
plored to great detail, and it has been investigated how 
and where cooperativity originates. In general, coopera- 
tive switches are much more robust against noise f^ ^ . 

The paradigm model for a cooperative genetic switch 
is the system made up by the host bacterium Escherichia 
coli which is infected by the parasitic bacteriophage T4 
(A-phage). A injects its own DNA into the host cell 
where the A-DNA fuses with the bacterium DNA. Thus, 
A is able to abuse the host cell facilities to either re- 
main dormant and get reproduced along with the bac- 



terium (lysogeny), or fuse new A-phages by the help of 
the host cell's miniature chemical plants, the process of 
lysis. The latter eventually leads to an array of a large 
number of new A inside E. coli. Finally the host cell bursts 
and releases a swarm of new A's pi 0]. Which of the 
two paths, lysis or lysogeny, is followed, is determined 
by a genetic switch which, in turn, is triggered by mes- 
senger molecules which we call repressor R (the protago- 
nist lysogeny agent) and aggressor R, the antagonist lysis 
agent. 

The typical, overall number of messenger molecules 
within a biological cell like E.coli is relatively small, rang- 
ing from a few to some 100. On the one hand, the chem- 
ical processes involved in the synthesis and degradation 
of the messenger molecules are noisy and, to some ap- 
proximation, described by a master equation [p|. Noise 
means that at some given time there is a surplus in the 
molecule synthesis in respect to the degradation, and the 
overall number of the respective molecule increases; and 
vice versa. That means that even in a state which is 
stationary on average, fluctuations in the number of the 
molecules occur. It has been extensively studied in how 
far the system can be influenced by such noise due to 
the feedback circle |, |, g, |, |, |l^, |ll|, [|. Note that 
this aspect of noise enters only in the number of involved 
molecules as a function of time. 

On the other hand, it might be argued that the spa- 
tial distribution of these molecules, governed by Brow- 
nian motion, may give additional cause for the influ- 
ence of noise. Usually, it is replied that the system is 
well-stirred, and therefore the molecules are always close 
enough to the switch on the DNA such that the spatial 
effects can be neglected. With the typical diffusion con- 
stant K r^ 2 ■ lO^^cm^/sec for a molecule of 50A size, 
the average diffusion time it takes to cross the cell is of 
the order of one to a few msec. Processes affected by 
genetic switching occur on time scales of the order of 
cell division of one or a number of cell generations, typ- 
ically tens of minutes or longer. The molecules, on this 



long time scale, are therefore well-mixed within the cell. 
The effects of noise on the switching process can, under 
this well-stirredness assumption, be considered directly 
in the biochemical circuit M, 0, M, pi Bl llG, , or it can 
be thought of as an activation process ||l3| . 

However, we will show that the involved system pa- 
rameters are not consistent with the above reasoning, 
and that well-stirredness alone is not sufficient to obtain 
a complete picture of the process. In contrast, we give ev- 
idence that the overall process consists of a large number 
of subprocesses during which there is an ongoing compe- 
tition between the protagonist and antagonist molecules 
which in turn gives rise to the influence of spatial fluc- 
tuations on genetic switching. Essentially, the reason for 
this claim is that the cell volume is large in comparison 
to both the size of the messenger molecules involved and 
the van der Waals (vdW) interaction radius around the 
operator sites on the DNA. If one divides the cell vol- 
ume into compartments of the size of this vdW radius, 
the occupation of individual compartments by the few 
molecules in the entire cell shows large fluctuations in 
time. 

In what follows, we exclusively consider the effects of 
spatial fluctuations in the above compartment picture. 
We distinguish the non-cooperative and the cooperative 
cases, and our model switch is supposed to work accord- 
ing to the following simple rules, (i) The non-cooperative 
switch gives rise to the state of lysogeny if one R is bound 
to the operator. This R can dissociate from the opera- 
tor with some time constant, and it can be replaced by 
another molecule, R or R, which is within the vdW inter- 
action volume (IV). If either another R substitutes the 
dissociated one, or the operator remains vacant, the dor- 
mant state is preserved. The switch is flipped, and lysis 
initiated, if eventually an R molecule binds to the oper- 
ator site, (ii) In the cooperative case, two R's can bind 
to the operator. In this configuration, the first facilitates 
the binding of the second. Only if both R's dissociate 
from the operator and are eventually replaced by one R, 
the switch flips towards lysis. Thus, after full dissociation 
of the R (or R's) from the operator, the question for both 
cases is whether there is at least one R and/or R within 
the IV. If only one species is present, we assume that 
binding of one molecule of this species necessarily occurs. 
If both species are present, we introduce a (1 — x) factor 
in favour of R binding (and x in favour of R binding). 

Consequently, our model can be rephrased as a re- 
newal process in the following sense. As the distribu- 
tion of messenger molecules outside the IV is irrelevant, 
the occupation of the IV can be regarded independent of 
the previous occupation after the diffusion time it takes 
a molecule to cross the IV, the renewal time St. With 
the typical vdW radius of lOOA, we find with the above 
K that St ~ 10^^ sec. Keeping track of the systems at 
"stroboscopic" times St, 2St, . . ., we can employ a simple 
statistical analysis. The basic ingredients are the prob- 



abilities Ho and Aq that neither R nor R is inside the 
IV, and that there is at least one such molecule present, 
1 — Ho and 1 — Aq. lip = {IV/cell volume} is the prob- 
ability that a single molecule is within IV, and there are 
-/Vr and A'p- molecules of either species within the cell, 
we find Ho = (1 - p)^« and Aq = (1 - p)^^ Q. 

Let us collect some relevant numbers. The radius of 
E.coli is of the order of l/xm, and the free volume in 
which the messenger molecules diffuse within the cell is 
^ l^iw? pi]. Comparing to the vdW radius, we obtain 
p ~ 5-10^^ Note that for these numbers, the probability 
that none out of 100 molecules is within the IV, (1 — 
p^ioo p^ 99.5%, is still very close to 1. The presence of 
such small numbers which give rise to the fact that the 
associated probabilities are either close to or to 1 is the 
reason for the relevance of spatial fluctuations. 

Consider first the non-cooperative case. Assume that 
the bound R molecule dissociates with the characteristic 
time scale r. Excluding that the molecule does not im- 
mediately rebind, one witnesses a competition between 
the two kinds of messenger molecules which can possi- 
bly bind to the relevant operator sites. This competi- 
tion is characterised through the four events Pi = HqAo, 

P2 = (i-no)Ao, P3 = (i-no)(i-Ao), P4 = no(i-Ao) 

which define the joint presence or absence of the two 
types of molecules. These four configurations can be sub- 
divided into those which leave the genetic switch in the 
dormant mode, i.e., which prevent an R molecule from 
binding to the operator site, and those which lead to R 
binding to its operator. The former comprise Pi and P2, 
the latter is given through P4. In turn, P3 defines a mixed 
state whose mean outcome will be (1 — x) in favour of R 
binding, and x in favour of R binding. The probability 
for inhibition is thus Pinhib = Pi + P2 + (1 — x)^3i the 
one for lysis is Piys = 1 - Pjnhib- 

The whole process can therefore be stripped down to 
the occurrence of a number i of inhibition events, termi- 
nated by a step leading to lysis, i.e., lysis will eventu- 
ally occur according to a sequence Pjnhib, Pinhib, ■ ■ • , Piys 
with joint probability Pinhib^iys ^^'^ normalisation A/" = 
Piys/(1 — Pinhib) — 1- The mean time to obtain lysis 
according to this diffusion picture is (St) — St X^^g {i + 
l)Piysi^*^hib' resulting in (St) = St/P\ys. The quantity 
(St) increases with growing A^r, with decreasing N^, or 
with decreasing x, as it should. (St) is the time due to 
the diffusion renewal process. To obtain the overall char- 
acteristic time for lysis, we have to add the binding times 
of order r. This delay can be included through the aver- 
age number of steps (i) = J^'zlo ^-Piys-Finhib = -Pinhib/Piys, 
weighted by the probability (P2 -I- (1 — x)^3) that a re- 
newal step actually involves a rebinding of an R. Multi- 
plied by T and added to (St) , this leads to the character- 
istic lysis time 
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which will be discussed in comparison to the time scale 
in the cooperative case. 

As mentioned, the cooperative scenario involves two R 
molecules. If one is already bound to an operator site, it 
facilitates the binding of another R molecule to the sec- 
ond operator site reserved for R such that % « 0. Usually, 
two R's are bound. The antagonist molecule R can only 
bind and initiate the divergence to the lytic track if both 
R sites are vacated. I.e., if the one R dissociates and does 
not rebind during the dissociation time of the second R. 
Moreover, one has to consider that not each time both 
R's are dissociated, R binds. In fact, some R molecule 
can bind to the operator sites and restart the dissociation 
process. As usually more R than R are within the cell, 
this case occurs more often, on average. Thus, if the char- 
acteristic time r^^ for the dissociation of the second R is 
large in comparison to the renewal time 5t, a sufficiently 
high number of R molecules makes it rather improba- 
ble that the R-related operator sites remain unoccupied 
long enough as to allow for the complete dissociation of 
R to occur: the characteristic time for lysis in the coop- 
erative case should be considerably higher than for the 
non-cooperative case Q . 

To quantify this cooperative process, let us assume 
that s — T^^ /5t is the number of renewal steps corre- 
sponding to the dissociation time of the second R that is 
still bound. After dissociation, 1 — Hq defines the prob- 
ability that, in one given renewal step, an R molecule 
binds to the vacant operator site and reconstitutes the 
initial configuration with two R's bound to the DNA. 
The probability that during 5t no such reconstitution oc- 
curs is given by IIq. The probability that reconstitution 
occurs in less than s renewal steps is then described by 
the combined process 77= (1 — IIq ) X)i=o -^0, obtaining 
rj = 1 — ng. The target process for finding the possibility 
for lysis thus corresponds to one of the following cas- 
cades of events, ?y, ????, . . ., ??'??, . . . where rj = 1 — rj. I.e., a 
certain number of superprocesses rj occurs during which 
reconstitution takes place, and finally no R replaces the 
dissociated first R until the second R dissociates, too. 
The associated mean number of "superprocesses" rj is 
(i)rf = v/{^ "V)- I^^ order to estimate the characteristic 
time connected to this process, we have to include two 
contributions. The first is the average time consumed 
by an rj superprocess; that is, r^ = 6t J2i=o (* + ^)^h — 
St{l - ng -h sng+^ - sng)/(l - Hq). The second is the 
dissociation time r^ elapsing after each re-binding of R. 
Finally, as not each complete dissociation of the two R 
molecules leads to a successful binding of the antagonist 
R molecule, we obtain the characteristic time scale 



j-lys 



-'lys-'inhib 



+ 



6t 

Piys 



(2) 



for the occurrence of cooperative lysis. In Eq. (0), the 
time constant for complete dissociation of both R's is 
given by {r^tl^) = {^)rrrrf+i{^)rf+^)T^- It is due to the 



additional weighting through {i)jf that T^^^ exceeds T^^^ 
considerably. 

Both characteristic times can be evaluated numerically. 
Essentially, the non-cooperative lysis time T^^, for a 
fixed number N'^, grows almost linearly in iVR, compare 
Fig. y. In contrast, T^^^ grows almost exponentially for 
fixed N^, eventually reaching extremely large values for 
higher iVR. For N-r, — 1, both characteristic times co- 
incide, as they should. Conversely, for fixed N-^, both 
characteristic lysis times fall off like a power law for in- 
creasing N^. 

These results are in qualitative agreement to those ob- 
tained from models considering exclusively number fluc- 
tuation: cooperativity enhances the accuracy of the sys- 
tem (the resistance against noise) exponentially, com- 
pare, e.g., Q. As our diffusion based model can lead to 
significantly large lysis times which can be of the same 
order of magnitude as the results from the number fluc- 
tuation models or even larger, depending on the assumed 
parameters, it may not be sufficient to consider number 
fluctuations only. In particular, it has been demonstrated 
that the well-stirredness assumption is no sufficient a pri- 
ori condition to exclude the spatial inhomogeneities aris- 
ing from the spatial diffusion of the molecules. This is 
based on the fact that each dissociation-rcbinding pro- 
cess is influenced by the fairly high probability that no 
molecules are in the interaction volume during a renewal 
step. 

It should be emphasised that our results are sensitive 
to the very numbers which are assumed for obtaining esti- 
mates for the characteristic lysis times. A small variation 
of these numbers can lead to a large change in the flnal 
result, so for a given system the parameters should be 
carefully verifled before estimates like the ones obtained 
herein are calculated. 

The basic ingredient of our model is the separation 
of the entire free cell volume into a bath constituted by 
the free diffusion volume, and into the IV. Due to this 
assumption, the very conflguration outside the IV can 
be neglected. Consequently, for Monte Carlo simulations 
of the combined process in which both the number of 
molecules and their spatial variation are random, the con- 
cept of the IV versus the free volume might prove useful 
in stripping off the unnecessary details and lowering the 
computation time. 

Our renewal time scenario relies on the existence of a 
more or less homogeneous distribution of the messenger 
molecules throughout the free diffusion volume such that 
the net exchange with the IV is approximately station- 
ary. In prokaryotic cells, this assumption should always 
be realised. It should also be valid in eukaryotes which 
feature a highly structured cell volume as long as there 
are no adsorption processes at cellular membranes which 
lead to immobilisation of the molecules according to a 
broad waiting time distribution which would give rise to 
a diverging exchange rate UM . 



By and large, in biophysics and biochemistry the role 
of noise in genetic circuitry, and cellular systems as a 
whole, has been increasingly assessed. This Letter shows 
that the spatial aspect of such fluctuations should not 
be neglected a priori, and its relevance for the process 
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should be checked. 
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we hnally remark that genetic switches are paradig- 
matic systems at hand of which effects of noise are stud- 
ied. The developed renewal-diffusion scenario therefore 
essentially pertains to a large variety of systems and pro- 
cesses involving lowly populated species which are spa- 
tially distributed and trigger foUowup processes on enter- 
ing some interaction zone, ranging from cellular feedback 
circles to clustering of bacteria, or to animal populations. 
In conclusion, there is no a priori well-stirredness condi- 
tion for such systems. 
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Here we assume that both kinds of molecules are of com- 
parable size. This is reasonable as both have to bind to 
similar operator sites. It is, however, straightforward to 
include differences in the size. 

Taking into consideration that a portion of the entire cell 
volume is occupied by other functional units like vesicles 
etc. 
Lysis can be prompted, for instance, by exposure of the 




Figure 1: Characteristic lysis times T^^^ and T^^'^ as function 
of A'^R. TJ'^ grows exponentially with Nr, i.e., it corresponds 
to an almost linear increase in the logarithmic case. Con- 
versely, T^Z'^ grows linearly. The drastic difference between 
these two patterns is demonstrated in the inset with linear 
axes. Parameters: 5t — 10~^sec, binding/dissociation times 
10"^sec, 7% = 5. 
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